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Abstract 

Wakes created by a parton moving through a static and infinitely extended quark-gluon plasma 
are considered. In contrast to former investigations collisions within the quark-gluon plasma are 
taken into account using a transport theoretical approach (Boltzmann equation) with a Bhatnagar- 
Gross-Krook (BGK) collision term. Within this model it is shown that the wake structure changes 
significantly compared to the collisionless case. 
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In ultrarelativistic heavy-ion experiments (SPS, RHIC) indications have been found for a 
temporary formation of a strongly coupled quark-gluon plasma (QGP) phase in the fireball 



created in nude— colons ,ee e,. Ref.fl). T h e main p robl e m of thi. p r o„ 
is to find clear signatures for the presence of a short-living QGP phase in this fireball. One 
promising class of signatures are hard probes, in particular high-energy partons of a few 
GeV or more, i.e. with an energy much higher than the temperature of the QGP. Energetic 
partons with a large transverse momentum, produced in initial hard parton collisions before 
the formation of an equilibrated fireball, propagate through the fireball. In this process 
they will lose energy, which depends on the state of the fireball (QPG or hadronic matter). 
Energetic partons manifest themselves at jets arriving in the detectors. The amount of jet 



' the jet through the fireball can 
3j. A strong jet quenching and 



quenching due to the propagation of the leading parton_o 
therefore serve as a signature for the QGP formation 
suppression of high hadron spectra have been observed at RHIC indicating the presence 
of a QGP phase 4J. Another phenomenon related to the propagation of a high-energy parton 
through the QGP are wakes and Mach cones which may be observable as conical flow and 
shock waves in particle correlations j^, [f], Q, [si, [9, Q. 

In electromagnetic plasmas wakes are a well known phenomena. Theoretically they can 
be investigated by considering transport equations. In the collisionless case they follow from 
the dielectric functions by solving the Vlasov equation together with the Poisson equation. 
It has been shown in this way that wakes are also created in a QGP, leading eventually to 
attractive potentials and Mach cones as well known in plasma physics (see below). Attractive 
potentials between partons might lead to bound states and may modify the J / ip suppression 
pattern 9|. In the static case, i.e. for a parton at rest, the wake potential reduces to the 
usual Yukawa potential describing Debye screening [l^ . 

Here we want to extend these investigations by taking into account collisions in the QGP. 
For this purpose we start from the Boltzmann equation approximating the collision ter m by 
the BGK description [3], allowing an analytic expression for the dielectric functions [l^.[l5j. 
In this approach the collision term is replaced by a momentum independent collision rate 
which we take as a parameter. This model has been used to study dispersion relations [ljj] 



and instabilities 



16| in a collisional QGP. Whereas the changes in the case of the dispersion 



relations are marginal - however the longitudinal plasma modes cross the light cone in a 
collisional plasma -, instabilities can be suppressed efficiently due to collisions. In the case 
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of complex pla smas, i.e., low-temperature discharge plasmas with microparticles such as 



dust grains [17], wake potentials due to ion flow, leading to an attractive 



brce between the 
]). Also collisions 



negatively charged microparticles, are investigated intensively (see e.g. 
between the ions and the neutral gas have been taken into account in this case and have 
been shown to truncate the wakes effectively 19[ . In general the collisions modify the wake 
potential by decelerating the ions. Therefore the positive ion cloud, responsible for the 
attractive part of the wake potential, is concentrated nearer to the dust grain, causing a 
narrower but deeper potential well. 

The longitudinal and transverse dielectric functions of a collisional QGP within the BGK 
approach are given by fl^ |. 
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where v is the collision rate and rri£> = \/l + Uf/6gT the Debye screening mass. Here ri/ 
is the number of quark flavors in the QGP. From the dielectric functions the dispersion 
relations follow 14j, e.g. the longitudinal mode (plasmon) from ei(uj,k) = 0. Note that uj 
has a real and imaginary part here. 

The induced charge density by a parton with color charge Q a propagating through the 



QGP with velocity v follows from [9| 
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Due to the delta function uj is restricted here to real values and is always in the space-like 
region for < v < 1. 



The induced charge density using cylindrical coordinates can be written as 201 ] 
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where J is the Bessel function, A = [Ree ; ] 2 + [Ime ; ] 2 , and k = \/k 2 — uj 2 /v 2 . 
The real and imaginary parts of the dielectric functions read 
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Solving together with (j3j) numerically [21] leads to the induced charge densities shown 
in Figs.l and 2. The induced densities are proportional to m 3 D and are scaled with mjj as 
well as with the color charge Q a of the moving parton. Assuming for example rif = 2, g = 2 
(corresponding to a s — 1/tt), and T = 250 MeV, we get a typical value of mjj = 580 MeV. 
In Fig.l the parton velocity is v = 0.55 c and v = 0.99 c in Fig. 2. The left panels show the 3 
dimensional plots for the choice v = to v = 0.8m£> 22j. The right panels show the contour 
plots of the induced charge density. One observes a clear modification of the induced charge 
density due to the collision rate. The general structure of the wakes, corresponding to a 
modified screening, a cone-like structure, and wave excitation has been discussed in detail 
in Ref.j^l in the collisionless case. However, in the case of collisions these structures are 
smeared out and less pronounced for increasing v. We interpret this as the fact that the 
collisions, leading to thermalization, reduce the anisotropy caused by the perturbation of 
the moving charge in the QGP. 

The wake potential is given by (see appendix) 

Qa i> ^—ik-(f—vt) 
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It is of the same form as in the collisionless case, where it can be shown that in cylindrical 



coordinates it reduces to 
$ a (p,z,t) 
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In Figs. 3 and 4 the wake potentials scaled with respect to mo are shown for a parton 
with a given color charge Q a for velocities v = 0.55 c and v = 0.99 c and various values 



for the collision rate from v — to v — 0.8rri£). It can be seen that collisions modify the 
strength and the structure of the wake potential significantly. Again, as in the case of the 
induced charge density the potential structure is washed out due to the collisions between 
the partons in the QGP. However, the potential becomes also stronger at finite collision 
rates, indicating that collisions in the QGP can enhance the attractive and repulsive parts 
of the interaction potential of a fast parton as in the case of non-relativistic plasmas. In 
particular, the attractive part behind the particle becomes deeper which might lead to an 
even stronger binding of diquark states propagating through the QGP as discussed in Ref. jo]. 
(Note that the induced charge density and wake potential diverge for point-like particles at 
the origin as in the static case corresponding to a Yukawa potential.) 

In conclusion, we have applied the BGK approach in order to study the influence of 
collisions on the wake formation in the QGP. Clearly, this approach allows only qualita- 
tive statements. A consistent QCD approach, in which for example the collision rate are 
momentum dependent, is beyond the scope of this work. Keeping this in mind, we regard 
the (momentum independent) collision rate in our model as a free parameter. Anyway, our 
model demonstrates clearly that collisions in the QGP lead to a significant modification of 
wakes created by fast partons compared to the collisionless case. In particular, the wake 
structure is less pronounced. Hence possible observable effects such as the conical flow and 
wave excitations will be reduced. On the other hand, the attractive potential well becomes 
deeper. Similar effects of collisions have been observed in electromagnetic plasmas. 
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I. APPENDIX 

Here we will explicitly derive equation ([6]) for the wake potential in the case of collisions. 
For the collisionless case it follows from a combination of the linearized Vlasov equation and 
the Poisson equation. Here we start from the kinetic equation with a BGK-collision term. 
The derivation from the equilibrium distribution function fo(p) is then given by (see (14) in 
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Ref.Q) 
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where q is the charge of the plasma particles and D := uj + iv — k ■ v with v = p/p in the 
ultrarelativistic case (v — 1). 
The quantity 77 is given by 



(9) 



N J r p dp 

with the equilibrium particle number density N and dp = d 3 p/(2-ir) 3 . 

The Poisson equation for the potential $ of a point charge Q with velocity v, where we 
drop the color index a for convenience, in the plasma reads 



A$(f, t] v) = -4nQ S(f- vt) - p(r, t), 



(10) 



where p — q J dp 8 f is the induced charge density in the plasma. In momentum space it 
reads 

k 2 $(LU,k,p) = 8n 2 Q5(LU-v- k) + q [ dp5f(u,k,'p). (11) 



Using E = — z/c$ in momentum space and dfo(p)/dp = vdfo(p)/dp we can write the Poisson 
equation by combining (TTTT) with ([8]) and Q) as 



k 2 $(uj, k, p) = 8tt 2 Q 8{u - v ■ k) + k 2 (l - X )$, 
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Solving (fT2"l) for $ and Fourier transforming it we obtain the wake potential 
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Comparing y with the longitudinal dielectric function in the BGK approximation (see (17) 
and (18) in [Jj), 

,2 



= 1 + JL.J dp{s . g).^. _ J| / dp(v -. J) /oD - 



(15) 



we find x( a; = v ■ k) = ei(uj = v ■ k) . Hence ( Fl4|) coincides with ([6]) and has the same form 



as in the collisionless case 



9]. 



This result can also be derived more generally, i.e. in an homogeneous and isotropic 
plasma with and without collisions, assuming linear response theory for the induced charge 
density 

j Pind 

(1/ei - l)pexu where the total charge density p tot = p ext + p ind , and 
using the Poisson equation $ = 4irp ext /(k 2 €i) with p ext = 2ttQ5(uj — v ■ k) after Fourier 
transformation. 
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FIG. 1: Left panel: Spatial distribution of the induced charge density for a parton moving with 
velocity v = .55c with collisions for (y = 0, 0.2m£>, 0.5m/), 0.8m/)) (Top to Bottom). Right panel: 
These plots show the corresponding contour plots 
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FIG. 2: Left panel: Spatial distribution of the induced charge density for a parton moving with 
velocity v = .99c with collisions for (y = 0, 0.2m£>, Oimu, 0.8m/)) (Top to Bottom). Right panel: 
These plots show the corresponding contour plots 
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FIG. 3: Left panel: Spatial distribution of the scaled wake potentials for a parton moving with 
velocity v = .55c with collisions for (y = 0, 0.2m/), 0.5m£>, 0.8m/) (top to bottom panel). Right 
panel: These plots show the corresponding equipotential lines. 
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FIG. 4: Left panel: Spatial distribution of the scaled wake potentials for a parton moving with 
velocity v = .99c with collisions for (y = 0, 0.2mo, O.brriD, 0.8mo (top to bottom panel). Right 
panel: These plots show the corresponding equipotential lines. 
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